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ABSTRACT 

We report the discovery of polarized flux at 2.2 /im from the bright shell of the w 320 year old supernova 
remnant Cas A. The fractional polarizations are comparable at 6 cm and 2.2 fim , and the polarization 
angles are similar, demonstrating that synchrotron radiation from the same relativistic plasma is being 
observed at these widely separated wavebands. The relativistic electrons radiating at 2.2 fim have an 
energy of « 150 GeV, (7 « 3 x 10 5 ), assuming an « 500 fiG magnetic field. The total intensity at 2.2 
/xm lies close to the power law extrapolation from radio frequencies, showing that relativistic particle 
acceleration is likely an ongoing process; the infrared emitting electrons were accelerated no longer than 
« 80 years ago. There is a small but significant concave curvature to the spectrum, as expected if 
the accelerating shocks have been modified by the back pressure of the cosmic rays; given calibration 
uncertainties, this conclusion must be considered tentative at present. The 2.2 /im polarization angles 
and the emission-line filaments observed by HST are both offset from the local radial direction by 10 
— 20 degrees, providing evidence that the magnetic fields in Cas A are generated by Rayleigh- Taylor 
instabilities in the decelerating ejecta. 

Subject headings: acceleration of particles — magnetic fields — polarization — radiation mechanisms: 
nonthermal — supernova remnants 



1. INTRODUCTION 

The energy distribution of « 0.1 — 10 GeV electrons in 
supernova remnants has been studied for decades through 
their radio synchrotron radiation. The distributions are 
typically power laws, with low energy spectral flattenings 
often seen due to absorption in the interstellar medium 
(Kassim 1989), or perhaps even due to internal ionized 
gas in the case of Cas A (Kassim et al. 1995). The power 
law slopes (spectral indices) vary significantly from rem- 
nant to remnant (Green 2001). Although it has long been 
accepted that diffusive shock acceleration is responsible for 
the GeV electrons (Bell 1978), there are still no clear ob- 
servational signatures for how this is regulated in different 
remnants. Spectral index variations within an individual 
remnant are often very small; this presents a dilemma be- 
cause the asymmetries in remnant structures and dynam- 
ics should lead to large variations in spectral index if par- 
ticle acceleration is due to first-order Fermi acceleration 
at shocks. In older remnants, such as the Cygnus Loop, 
spectral variations probably reflect local differences in ab- 
sorption and confusion along the line of sight (Leahy & 
Roger 1998). But in Cas A, the spectral variations across 
the remnant (Anderson & Rudnick 1996) are most consis- 
tent with different power laws, indicating local variations 
in relativistic particle acceleration (Wright et al. 1999). 
Information about deviations from the power law shape 
could provide important information about the acceler- 
ation mechanisms and the possible modifications of the 
shocks by cosmic-ray backpressures (Berezhko & Ellison 
1999); only weak evidence for such deviations exists to 



date (Reynolds & Ellison 1992). 

The maximum energy to which cosmic rays can be ac- 
celerated in supernova remnants is another issue of cur- 
rent importance. Non-thermal X-ray synchrotron radia- 
tion from the supernova remnants SN1006 (Koyama et al. 
1995; Reynolds 1996), RX J1713.7-3946 (Koyama et al. 
1997; Slane et al. 1999), Cas A (Allen et al. 1997), and 
RCW86 (Borkowski et al. 2001) suggest that relativistic 
electrons have energies up to 100 TeV. Gamma-ray detec- 
tions of SN1006 (Tanimori et al. 1998) and RX J1713.7- 
3946 (Muraishi et al. 2000) have been interpreted as the 
Inverse Compton emission from these high energy elec- 
trons against the cosmic microwave background. 

Although the energies of electrons that would produce 
synchrotron radiation in the infrared are considerably 
lower than those emitting X-rays and 7— rays, they offer 
the possibility of a model-independent probe of electrons 
at much higher energies than seen in the radio. The pos- 
sibility of near or mid-IR synchrotron emission was first 
suggested to us by R. Tuffs (private communication and 
Tuffs et al. (1997)), based on their impressions of the 
ISOCAM images (Lagage et al. 1996). The possibility was 
again raised by Gerardy & Fesen (2001), who noted the 
close similarity between their 2.2 /*m (but not 1.2 /jm) im- 
age and the radio images. This led us to seek verification of 
the infrared synchrotron origin through polarization mea- 
surements. 

2. OBSERVATIONS 

2.1. Choice of Field 
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The observations were made using NSFCAM on the 
IRTF in polarimetry mode on 1 August, 2001. Since NSF- 
CAM has a field of view of about 75" (0.3" per pixel) and 
we would have time to observe only one field, we had to 
make a careful choice of location within the Cas A rem- 
nant. The 2.2 //m band image kindly provided by C. Ger- 
ardy and R. Fesen (Gerardy & Fesen 2001) gave us infor- 
mation on the flux from Cas A at 2.2 /zm . Our own radio 
images of Cas A (Anderson & Rudnick 1996; Anderson, 
Keohane & Rudnick 1995; Gotthelf et al. 2001) provided 
us with total intensity and polarized intensity maps that 
could be compared to the 2.2 jim band image. 

Our choice of which region in Cas A to observe was in- 
fluenced by a number of factors. First, there had to be a 
good correspondence between the radio synchrotron emis- 
sion and the 2.2 /im emission. This was done to avoid small 
bright regions of emission at 2.2 /im that might be contam- 
inated by emission lines unrelated to the radio synchrotron 
emission. Second, there had to be significant polarized in- 
tensity in the radio, which we hypothesized would corre- 
spond to similar polarization in the near infrared. Third, 
there had to be sufficient flux at 2.2 /im to make fractional 
polarization measurements of a few percent, a level similar 
to the observed radio polarization. Using these criteria, we 
chose a segment of Cas A's bright ring in the northwest 
(azimuths of w Oto — 40 degrees) centered at RA 23h 23m 
26s , Dec 58°, 48' 58.4", (J2000), as shown in Figure 1. 

2.2. Polarization Observations 

The techniques for using NSFCAM as an imaging po- 
larimeter are described in Jones (1997, 2000). Images are 
taken of the source and an off-source sky location at four 
positions of a half waveplate that can be rotated in the 
beam. This allows us to form sky subtracted intensity im- 
ages at polarization position angles of 0, 45, 90 and 135° 
(waveplate angles of 0, 22.5, 45 and 67.5°). We formed two 
independent estimates of the total intensity (Iq + Igo) & n d 
(/45+/135), and found that these were consistent with each 
other within the off-source noise errors. Photometric cal- 
ibration was done using observations of HD 225023 from 
the Elias standards list (Elias et al. 1982). The 2.2 /im (K) 
intensity scale was set using the conversion 0.0 uik — 635 
Jy. The final accuracy of our spectral indices is limited by 
the 2.2 /im extinction corrections, as discussed below. 

From the images at four angles, we then computed 
Stokes images Q (Iq — /go) and U (-Z45 — /135) and the 
total intensity image |(Io + I45 + ho + ^135)- Note that 
these Stokes parameter definitions differ by a factor of two 
from those commonly used in the radio. The details of our 
sequencing of source, sky and waveplate position differed 
somewhat from Jones (1997). For Cas A we rotated back 
and forth between two waveplate positions (either for Q or 
U) for several minutes, then repeated the same procedure 
on the sky position. Integrations at a single waveplate po- 
sition were 10 seconds. It took about 2 seconds for the 
waveplate to move to a new position. 

NSFCAM is not a true imaging polarimeter, but rather 
a camera retrofitted with some polarization optics, so we 
were limited by systematics in the precision of our po- 
larimetry to about ±0.3%. This limit is due primarily to 
fluctuations in atmospheric transmission between the time 
it takes to integrate at different positions of the waveplate. 



Since we are interested in polarizations of a few percent 
or more in Cas A, this instrumental limitation did not 
hinder our observations. In fact, the nebulosity in Cas A 
at 2.2 /im is sufficiently faint that our measurements are 
effectively photon noise limited by the background. 

Calibration of the position angle was done using obser- 
vations of AFGL 490 two nights earlier with the same ob- 
serving procedure as for Cas A. AFGL 490 was assumed to 
have a polarization angle of 9 — 115° at 2.2 /im (Kobayashi 
et al. 1978). The Q' and U' values quoted in Table 1 are 
given as observed; the PA contains a —33 degree correction 
to place it in the proper sky frame. An extra check on the 
position angle calibration was done using observations on 1 
August of the reflection nebulosity surrounding the young 
star RNOl (Weintraub & Kastner 1993). The centro- 
symmetric scattering pattern around this star produces 
all possible position angles and provides a very sensitive 
calibration. Calibration of the efficiency (which exceeds 
98% at 2.2 /im ) was performed in previous observing runs 
using SI in p Oph which has a polarization of 1.9% at 2.2 
/im . Instrumental polarization was checked during the 
night by observing an unpolarized standard star from the 
UKIRT faint standards list (FS 149) and was found to be 
too small to be measured. 

The telescope was offset 100" to the North for the sky 
frames. This was intended to move the telescope well off 
Cas A, but not so far that telescope movement and set- 
tling time seriously reduced our observing efficiency. In an 
effort to increase our time on Cas A itself, we reduced the 
number of background frames to about half the number 
of frames taken on the source. The total on source inte- 
gration time was ~ 1 hr. Serious telescope drift during 
the long integrations on Cas A resulted in the movement 
of stars in the background frames relative to the source 
frame. This produced elongated stars in the background 
subtracted images (see Figure 1,2.2 /im intensity). No cor- 
rection has been made for these problems, and the effect 
on our spectral measurements is described below. Our po- 
larization measurements are less affected by the star trails 
than the total intensity because the star polarizations are 
expected to be « 1% as discussed below, while Cas A's 
polarization is w 5 — 10%. 

Since the time between source and background frames 
ranged over several minutes, fluctuations in sky back- 
ground could cause a residual DC level to be present in 
the sky subtracted images. Although small compared to 
the net intensity of the brighter nebulosity in Cas A, this 
DC offset could be comparable to the polarized intensity, 
and produce a spurious net intensity in the Q and U im- 
ages. These offsets were removed by first identifying sev- 
eral small regions of the sky-subtracted total intensity im- 
age with values fluctuating near zero for the intensity. We 
simply assume the true intensity is zero when averaged 
over these locations and subtract the respective means 
from the I, Q and U images to force them to be zero there. 
This procedure will remove the DC offset produced by sky 
background fluctuations in the Q and U images, but also 
would remove any real, but faint emission. Sensitive ob- 
servations in the radio do show faint total intensity and 
polarized emission from these regions, at intensities of ap- 
proximately 15% that of the bright ring; no correction has 
been made for this small bias. 
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2.3. Radio and Optical Images 

The radio images were produced from observations at 
the Very Large Array 2 using all four configurations over 
the time period from May, 1997 through March, 1998. 
Maps were made by using the data from 1.285 GHz for 
the 20 cm band, and from 4.64 GHz for the 6 cm band. 
The flux calibration scale was set using standard VLA pro- 
cedures (VLA Calibration Manual 2002) and the source 
3C48. The polarization calibration also used standard 
VLA procedures with the angle fixed by assuming a po- 
sition angle of —10° for 3C138. Other details of analysis 
and image creation are similar to those described by An- 
derson, Keohane & Rudnick (1995). The full resolution of 
these images was w 1.5". 

The optical emission in Cas A's bright ring is domi- 
nated by clumpy emission lines on a wide range of scales, 
most recently studied in great detail by Fesen et al. (2001). 
The emission line material represents a distinct tempera- 
ture component, separate from the synchrotron plasma, 
but which may dominate the dynamics driving the rela- 
tivistic particles and magnetic fields. In particular, Fesen 
et al. (2001) noted from their high resolution WFPC2 ob- 
servations, a very large number of filamentary structures 
that they suggested were Rayleigh- Taylor instabilities. Al- 
though there is significant small-scale scatter, we measured 
the approximate position angle of the filaments at the az- 
imuths of the boxes used for the synchrotron measure- 
ments and list these in Table 1. 

3. DATA ANALYSIS & RESULTS 

3.1. Image Registration 

The registration of the 2.2 /zm and radio images was 
done in two steps. The radio coordinates from the VLA 
are already in the FK4 reference frame to an accuracy 
of s=s 0.1". To put the 2.2 /zm image in the same frame, 
we first used the Sloan Digital Sky Survey to register the 
image of the entire Cas A field kindly provided by C. Ger- 
ardy & R. Fesen (Gerardy & Fesen 2001). We then mini- 
mized the rms differences between our new small polarime- 
try field and the Gerardy/Fesen image, allowing for shifts 
in position, rotation, and scale. The final accuracy of this 
registration is < 0.1", and the uncertainties have no effect 
on the results discussed here. 

3.2. Spatial averaging 

The signahnoise on the 0.3" pixel level was not sufficient 
either for accurate spectral measurements or determina- 
tion of the polarization. We therefore found the average 
brightness for I, Q, and U in each waveband in six boxes 
as shown in Figure 1, and listed in Table 1. These boxes 
were chosen to have sufficiently high polarized flux at 6 
cm, to avoid regions with bright stars in the sky ("off") 
frame, and to avoid the edges of the frame where the flat 
fielding is less accurate. 

For display purposes, we also convolved the polariza- 
tions (Q and U) to a resolution of 10". The resulting low 
resolution polarized intensities arc shown on the right side 
of Figure 1. The middle image of the right side represents 
a prediction of the 2.2 /zm polarizations using the radio 

2 The Very Large Array is an instrument of the National Radio Astrc 
cooperative agreement by Associated Universities, Inc. 



data and 2.2 /im total intensities. It was constructed by 
taking the fractional polarization image at 6 cm and mul- 
tiplying it by the total intensity image at 2.2 /zm . The 
overall correspondence is good, but not perfect. The 2.2 
/zm polarized intensity map contains, e.g., a bright spot 
due to residuals from the bright star immediately to the 
SE of box 6. In addition, it has not been corrected for the 
noise bias, which is responsible for some of the structure 
in the image. The 6 cm polarization angles were used, af- 
ter correction for Faraday rotation as described below, as 
a prediction for the 2.2 /zm angles. The correspondence 
here is somewhat better than expected from the errors. 

A number of corrections to the total and polarized in- 
tensities are needed before the 2.2 /im and radio intensities 
and polarizations can be quantitatively compared. First 
is the contamination of the 2.2 /zm total intensity by fore- 
ground stars in the boxes. The fluxes of the brightest 
stars in each box were therefore subtracted yielding the 
"corrected" 2.2 /im intensities given in Table 1. We es- 
timate that the remaining positive star contamination is 
less than a few percent. 

More troublesome are the negative trails from stars in 
the sky frame, since the contribution from each individ- 
ual star is spread across the box, and is therefore fainter 
and harder to recognize. Although the boxes avoid the 
brightest negative trails, there are certainly unrecognized 
trails that reduce the 2.2 /im intensities. In individual re- 
gions, such as the NW corner of the frame, we can look 
for negative trails in the absence of bright diffuse emis- 
sion from Cas A. The distribution of pixel values in that 
region is consistent with a gaussian distribution around 
mean zero with a long negative tail. For this region, the 
negative tail can be recognized and ignored in determin- 
ing the zero offset level. When negative star trails overlap 
bright 2.2 /zm emission, however, they can reduce the in- 
tensity without being recognized. There is thus a negative 
bias in the 2.2 /im intensities, which we expect to be small 
on average (« 0.5/zJy/D", or « 5% of the total intensity), 
but could possibly be as much as s=s 30% in an individual 
box, without being recognized. The 2.2 /zm total inten- 
sity values quoted in Tabic 1 could therefore actually be 
somewhat larger, but because of the scientific significance 
of such larger values, we have taken a more conservative 
approach and not tried to make any corrections for this 
possible bias. 

Another bias correction that we did not make to the data 
is due to the fact that the sensitive radio observations show 
that synchrotron emission is present everywhere in the 2.2 
/zm image, and therefore it is not really possible to define 
an appropriate zero, as discussed earlier. The characteris- 
tic intensity of the radio plateau is approximately 15% of 
the bright ring, so defining it to be zero at 2.2 /zm results 
in spectral indices between 2.2 /zm and 6 cm that are too 
steep by « 0.01. 

Cas A's high optical extinction (> 4 m Searle (1971) ) 
requires that corrections are needed even at 2.2 /im . Hur- 
ford & Fesen (1996) find reddenings of 4.6 < Ay < 5.4 
magnitudes for five FMKs and Ay — 5.3 and 6.2 magni- 
tudes for two QSFs, and provide an extensive discussion 
of the literature on Cas A's extinction. In the west, Ay 

' Observatory, a National Science Foundation facility operated under 
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may extend up to 8 magnitudes, according to the CO and 
HI analysis of Troland, Crutcher & Heiles (1985). In this 
paper, we use a value of Ay = 5 magnitudes, which leads 
to A2.2fj.rn = 0.55 magnitudes, adopting the conversion fac- 
tors of Bessell, Castelli, & Plez (1998). Application of this 
extinction value results in the lines labeled "de-reddened" 
for I, Q, and U in Table 1, and the values plotted in Figure 
2. We also show, in Figure 2, the range of corrections that 
would result for the range 4.6 < Ay < 6.2 magnitudes. 

For any reasonable value of the extinction, it thus ap- 
pears that the 2.2 /im data consistently fall above the ex- 
trapolation from radio wavelengths, giving the spectra a 
concave shape, i.e., flatter at short wavelengths. There 
are several caveats to this conclusion. First, a mixture of 
different power-laws, as observed at cm wavelengths (An- 
derson & Rudnick 1996; Wright et al. 1999), will result 
in a concave spectrum. The rms variations within our 
boxes arc <= 0.01, insufficient to cause a significant ef- 
fect at 2.2 /im . However, it is likely that there is some 
contamination of the bright ring with emission from the 
steeper spectrum plateau, at least at radio wavelengths. 
The maximum contamination, for plateau emission with 
a icm™ = ^ 0.9 would result in an increase of the 2.2 /im 
flux density of « 35%. The actual 2.2 fim contamination 
is likely to be considerably lower, given our procedure for 
setting a zero background, as discussed earlier. However, 
in any case, this possible contamination is much less than 
the factor of two enhancement seen at 2.2 /im after cor- 
rection for extinction. 

A second concern is our calibration to the radio flux 
scale (VLA Calibration Manual 2002), which itself is 
based on "absolute" calibrations of Cas A (Baars et al. 
1977). These two steps are expected to be accurate to 
« 2%, and so could each contribute another 25% error to 
the 2.2 /im extrapolation. We therefore consider our detec- 
tion of spectral concavity to be reasonable, but somewhat 
tentative until the extinction and flux scale issues can be 
fully resolved. At present, we are analyzing the 6 cm — 2.2 
/zm spectra across the whole remnant, which should allow 
us to use differential measurements to separate real from 
calibration effects. 

3.3. Polarization 

The 2.2 /zm Stokes parameter values in Table 1 are listed 
as Q', U'; these represent the actual observations, which 
are only nominally Q and U prior to the calibration of 
the absolute polarization angle. The properly calibrated 
polarization position angles are listed in Table 1. 

In either the (Q',U') or (Q,U) frame, the observed polar- 
ized intensity {\f[Q 2 + U 2 )) represents a biased measure 
since it would give a positive result from noise, even in the 
absence of a true signal. As a simple correction for this 
bias, the fractional polarizations for 2.2 /tm listed in Table 
1 were calculated as %P = [y^iQ 2 -u'q+U 2 -of/-)]/ / where 
were calculated as the rms scatter among boxes of 
approximately 35 pixels on a side, away from bright stars 
and away from the bright ring. 

For the 2.2 /im polarization measurements, we also need 
to consider the possible contamination by the foreground 
interstellar polarization expected from the large extinc- 
tion towards Cas A. Catalogs of optical polarimetry (e.g. 
Heiles (2000)) have no stars in the same field as Cas A, 



so a direct measurement of the interstellar polarization in 
the visual has not been made. We have measured the po- 
larization of the bright star in the northwest of our image 
and find a 3a upper limit of 0.9% at 2.2 /im , although we 
do not know if it is as far away as Cas A. The expected 2.2 
/xm polarization for Cas A, given the visual extinction of 
5-6 magnitudes, can be estimated using the mean trend in 
polarization with extinction analyzed in Jones (1989) and 
Jones, Klebe & Dickey (1991) of' 

< P2.2pm(%) >= 2-2 * T 2 . 2fl m = 0.2 * A V RJ 1%. 

Thus, we conclude that interstellar polarization is not a 
significant contaminant of our polarimetry of the nebulos- 
ity at 2.2 /xm . After consideration of all these effects, 
and the faintness of the nebulosity being measured, the 
errors in the fractional polarizations and angles at 2.2 /im 
are dominated by the random contribution from sky back- 
ground photons, as opposed to systematic effects. The 
measured percentage polarizations at 2.2 /im range from 
« 4 — 10 %, comparable to those at 6 cm. 

To compare the 2.2 /im polarization angles to those at 
radio wavelengths, we must first correct for Faraday Ro- 
tation, which can arise both local to Cas A and along the 
line of sight. At 20 cm, the rotation is strong enough to 
actually depolarize the emission (Anderson, Kcohane & 
Rudnick 1995), and render rotation measure (RM) mea- 
surements at that wavelength unreliable for the higher fre- 
quencies. Kenney & Dent (1985) present RM measure- 
ments at low resolution but high frequencies, with a value 
of —106 rad mT 2 for the northwestern part of the rem- 
nant. In order to look for variations in RM among the 
boxes used here, we created three different polarization 
maps at 4.41, 4.64 and 4.99 GHz, using the VLA D con- 
figuration data alone to produce a beam size of w 10". 
The signaknoise was insufficient to detect significant varia- 
tions among the boxes, but averaging over all boxes yields 
RM — —105 ± 3 rad m~ 2 , consistent with the Kenney 
& Dent (1985) value. We therefore adopt the value of 
— 106 rad mT 2 for all boxes, and show the Faraday cor- 
rected position angles in Figure 1 and in the angle com- 
parisons at the bottom of Table 1. 

In Figure 3, we compare the position angles of the mag- 
netic fields derived from the polarization vectors in each 
box with the approximate angles of the filamentary struc- 
tures in the WFPC2 images of Fesen et al. (2001). All 
of the position angle measures change systematically with 
azimuth, with a small but significant clockwise offset from 
the local radial direction at least for the 2.2 /im and fila- 
ment angles. It is not clear if the 6 cm angles show any 
offset, but better local rotation measures would be needed 
to make the comparison with 2.2 /zm angles more accu- 
rately. 

4. SCIENTIFIC IMPLICATIONS 

4.1. Maximum energy of accelerated electrons 

Any mechanism to accelerate relativistic particles yields 
some maximum particle energy. The energy can be limited 
by a variety of factors including the time available for ac- 
celeration, the escape of the highest energy particles from 
the finite sized acceleration region, or the balance between 
energy gains and losses (Reynolds 1998, 2001). In Cas A, 
a hard X-ray tail, which falls two orders of magnitude be- 
low the flux extrapolated from radio wavelengths, suggests 
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that the synchrotron emission may extend to that regime 
(Allen et al. 1997). This conclusion has been called into 
question by Blcekcr et al. (2001), who argue that the 10 
— 15 keV emission is diffuse, and not concentrated at the 
outer shock as it should be if it were synchrotron. This 
argument is questionable, however, because the radio syn- 
chrotron radiation also does not show an enhancement at 
the location of the outer shock; Xray synchrotron emis- 
sion is therefore still a possibility. Bremsstrahlung pro- 
vides an alternative explanation for the hard X-ray emis- 
sion (Laming 2001), where lower hybrid waves in Cas A's 
strong magnetic field scatter electrons up to energies of 
tens of keV. There is not yet any good test to distinguish 
between bremsstrahlung and synchrotron explanations for 
Cas A, so the nature of the 10 - 100 keV emission is there- 
fore still open to question. With the demonstration of a 
synchrotron origin for the 2.2 /im emission, it is therefore 
interesting to estimate the maximum energy of the accel- 
erated particles. 

We start by calculating the value of the field strength 
that minimizes the energy in both relativistic particles and 
magnetic field (B m i n ). We assume a pathlength through 
the emitting regions of 0.5 pc, with a filling factor of unity. 
We assume no energy in relativistic protons. These as- 
sumptions lead to a field strength of w 500 /iG for the 
regions in Table 1. This field strength can then be used 
to estimate the energy of the relativistic electrons that are 
radiating at 2.2 /im . In a « 500 /iG field, this leads to 
energies of ~ 150 GeV, or a relativistic gamma of 3 x 10 5 . 
As argued below, the synchrotron spectrum almost cer- 
tainly extends a factor of 10 higher in frequency, which 
would result in particle energies(gammas) of « 450 GeV 
(9 x 10 5 ). This is still a factor of at least 10 3 below ener- 
gies to which SNRs are commonly expected to accelerate 
cosmic rays, but it is the first unambiguous demonstration 
of the extension of the same electron population to such 
high energies in a shell SNR. 

4.2. Synchrotron loss limits 

We use the same field strength calculation, and the lack 
of spectral steepening between 20 cm and 2.2 /im to place 
a limit on the time since the relativistic particles were last 
accelerated. The assumptions in the field strength cal- 
culation bias the lifetimes to longer values, so our lifetime 
estimates should be considered upper limits. At a factor of 
10 in frequency below the "break" frequency, the spectrum 
should already have steepened by 0.1 (e.g., Leahy (1991), 
Figure 3.3a). Since, if anything, the spectra have actually 
flattened at high frequencies, we adopt 0.1 as a conserva- 
tive upper limit to the steepening at 2.2 /im , yielding a 
"break" wavelength of < 0.22/im. This leads to lifetimes 
for the electrons radiating in the infrared that are no longer 
than 80 years. It is thus likely that particle acceleration 
is an ongoing process, at the current stage of Cas A's evo- 
lution. Based on our latest dynamical picture, this is the 
period when Cas A has probably swept up « 1 — 10 times 
its own mass (Gotthelf et al. 2001; Willingalc et al. 2002; 
DcLancy & Rudnick 2002) and during which the reverse 
shock is slowly coming to the end of its outward motion 
(in the sky frame). The reverse shock is still strong at 
this stage, and is encountering outward moving ejecta at 
speeds of « 5000 - 6000 km/s (Reed et al. 1995). The 



x-ray emitting filaments are rapidly decelerated to « 3500 
km/s while those in the radio are slowed to « 1800 km/s 
(Koralesky et al. 1998; Vink et al. 1998; DeLaney & 
Rudnick 2002); this provides at least one possible source 
of energy for the relativistic particle acceleration. 

4.3. Non-linear shocks 

If the acceleration of cosmic rays at SNR shocks is ef- 
ficient, then there is a back reaction of the cosmic rays 
on the shock strength and structure (Drury & Volk 1991). 
When the maximum cosmic ray momentum is limited by 
shock properties such as geometry, the compression ratio 
can become arbitrarily large (Eichler 1984). The increas- 
ing gyro-radii at higher energies allows those particles to 
see a higher compression due to the sub-shock structure, 
so the resulting spectrum becomes concave (Bell 1987; 
Berezhko & Ellison 1999; Reynolds & Ellison 1992), as ob- 
served here. The radio spectra can then be much steeper 
than the —0.5 expected in the strong shock test parti- 
cle limit, and can flatten out considerably before they cut 
off, e.g., in the X-ray, due to radiative losses and escape. 
Reynolds & Ellison (1992) suggested that slight curvature 
might have been observed in the integrated radio spec- 
tra of Tycho and Kepler SNRs, although the signal was 
marginal. 

In this work, we have shown that the 2.2 fim data, after 
correction for extinction, fall significantly above the power 
law extrapolation from the radio. The mean flattening be- 
tween o^ocm an d a 6crn, m 1S 0-06 (Table 1), equivalent to 
higher fluxes at 2.2 /im by a factor of « 2 . If the curva- 
ture were constant across this full range of wavelengths, 
then the local spectral index would change from w —0.75 
at cm wavelengths to ps —0.5 at 2.2 /im . This is essen- 
tially the same as the model shown by Berezhko & Ellison 
(1999). Using their Figure 5, we made a rough calculation 
of the flux expected at 2.2 /im for Box 1; this is shown as 
the nested square symbol in Figure 2. The agreement is 
quite good, suggesting that the signature of acceleration 
at modified shocks has been observed. In the context of 
the Berezhko & Ellison (1999) model, this implies a criti- 
cal injection rate of relativistic particles of >10~ 4 and, for 
most of the parameters they explored, efficiencies of 10 — 
100 %. Thus, with the caveats described earlier regarding 
the extrapolation of the radio spectrum, relativistic par- 
ticles appear to be an important part of the dynamics of 
the shock and the overall evolution of the SNR. 

4.4. Polarization 

The fractional polarizations and the polarization angles 
are very similar at 6 cm and 2.2 /im . This provides con- 
firming evidence not only for synchrotron radiation, but 
that we are looking at the same population of electrons 
at both wavelengths, validating the detailed spectral com- 
parisons above. The 2.2 /im results confirm the expec- 
tation from the Kenney & Dent (1985) mm wave studies 
that fractional polarizations of « 5 — 10% represent real 
field disorder, as opposed to Faraday depolarization. Thus, 
although we speak about the "radial" magnetic fields in 
Cas A, this only refers to the small net ordered component 
field. Approximately half of the field energy must be in a 
component which is disordered on scales < 1" (< 0.02 pc 
at a distance for Cas A of 3.4 kpc). 
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Gull (1975) first suggested that Rayleigh- Taylor insta- 
bilities in the decelerating bright ring of Cas A and other 
young SNRs would lead to radial magnetic fields. Detailed 
numerical simulations since that time confirm the preva- 
lence of these filamentary structures and their role in the 
amplification and alignment of the magnetic fields (Jun & 
Norman 1996; Jun & Jones 1999). Approximately radial 
filamentary structures have been seen in the Crab SNR, 
and, more recently in the WFPC2 images of Cas A (Fesen 
et al. 2001). In the northwest region investigated here and 
shown in Fesen et al. (2001) Figure 8, the Cas A filaments 
emerge from the suggested location of the reverse shock in 
directions that are systematically non-radial. These non- 
radial motions could arise from large scale asymmetries 
in the ejecta, as is well established from the distribution 
of different elements seen in the X-ray (e.g. Hwang, Holt 
& Petre (2000)). Alternatively, gradients in the circum- 
stellar medium could cause asymmetric deceleration (e.g., 
Koralesky et al. (1998)), again leading to large scale veloc- 
ity shears. In any case, the close correspondence between 
the R-T finger orientations and the local magnetic field di- 
rections strongly supports a connection between the two. 
Blondin & Ellison (2001) show that the R-T instabilities 
generated at the reverse shock can extend all the way out 
to the outer shock front, providing radial magnetic fields 
instead of the tangential ones otherwise expected at the 
shock. The transition region between the radial and tan- 
gential regions may have been observed in Cas A (Gotthelf 
et al. 2001), although the signal-to-noise is quite low. 

5. CONCLUSIONS 

The extension of Cassiopeia A's known synchrotron 
spectrum to 2.2 /zm opens up a variety of important issues 



regarding the acceleration of relativistic particles. First, 
we have established that the acceleration mechanism, usu- 
ally assumed to be diffusive acceleration at shocks, oper- 
ates to energies of at least 150 GeV in Cas A. In addition, 
there is reasonable evidence that these shocks are efficient 
particle accelerators, arising from non-linear interactions 
between the particles and the shocks. This creates the op- 
portunity to investigate these interactions in conjunction 
with the dynamical information about the several coupled 
plasmas (thermal and non-thermal) in the remnant. In 
turn, this will strongly constrain models for the high X- 
ray energy tail in Cas A's spectrum and assist in the chal- 
lenge of isolating the energy in relativistic protons. To 
make progress on these issues, a detailed investigation of 
the radio/2.2 /xm spectral shape as a function of the radio 
spectral index variations across the remnant is underway. 
The association of magnetic fields with the local R-T insta- 
bilities opens the door for investigating the accompanying 
amplification of magnetic fields in the remnant, and the 
roles of both relativistic particles and fields in the further 
dynamical evolution of the remnant. 

We appreciate the 2.2 /xm images provided by C. Ger- 
ardy and R. Fesen, and a number of important discus- 
sions on the science with Tom Jones. Support for work 
on Cas A at the University of Minnesota was provided 
by National Science Foundation grant AST-96- 19438 and 
for comparisons with HST data by NASA through grant 
number HST-AR-090537.01-A from the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. The comments of the referee were 
very helpful in improving the paper. 
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Table 1 



Intensity and Polarization Measurements 





I > . . i 

COX I 


Box I 


i > < i 

Box 3 


Box 4 


1 > r 

Box 5 


Box b 


Units (Error) 


Azimuth 


-10 


-20 


-30 


-35 


-42 


-32 


degrees 


hum observed 


12.0 


18.6 


7.6 




13.5 


12.8 


/iJy/D" (0.26) 


h^m corrected 


11.0 


9.6 


7.6 




9.3 


12.8 


/xJy /□" (0.26) 


hfim de-reddened 


18.3 


15.9 


12.6 




15.4 


21.2 


uJy/D" (0.43) 


Q'lam observed 


-0.67 


-0.48 


-0.54 


-0.84 


-0.84 


-0.48 


AtJy/D" (0.2) 


Q'lum de-reddened 


-1.11 


-0.80 


-0.90 


-1.4 


-1.4 


-0.8 


liJy/n" (0.3) 


U^m observed 

de-reddened 


-0.34 


.041 


0.33 


0.42 


0.50 


1.0 


AtJy/D" (0.2) 


-0.56 


.068 


0.55 


0.70 


0.83 


1.66 


/iJy/D" (0.3) 


% Pl^im 


6.2 


3.8 


7.2 




9.9 


8.3 


(1.5) 


PA^ 


70 


55 


41 


44 


42 


57 


degrees (13) 


hocm 


56,800 


53,800 


63,900 


56,300 


63,300 


61,600 


nJy/D" (300) 


hem 


21,500 


20,400 


24,200 


21,500 


23,500 


23,200 


AtJy/D" (100) 


Qficm 


-144 


51 


258 


780 


1535 


163 


MJy/LT 7 (45) 




1400 


965 


1550 


1270 


700 


680 


M Jy/D" (45) 


% Pecm 


6.5 


4.7 


6.5 


6.9 


7.2 


3.0 


(0.2) 


P A§ cm 


48 


43.5 


40.3 


29.2 


12.2 


38.2 


degrees (2) 


«PA (filaments) 


-40 


-40 


-33 


-70 


-78 


-33 


degrees (15) 


COMPARISONS 
















% Pl.1jiml%P%cm 


0.95 


0.81 


1.11 




1.38 


2.8 


(0.2) 


PA 6cm (corr) - PA 2fi 


-10 


-1.5 


+ 11 


-3 


-18 


-7 


degrees (13) 


6cm 
a 20cm 


-0.76 


-0.75 


-0.76 


-0.75 


-0.77 


-0.76 


(.01) 


2.2/im 
a 6cm 


-0.69 


-0.70 


-0.73 




-0.71 


-0.68 


(< -oi) 



Note: The absolute position angle (PA) on the sky at 2.2 ^im is (0.5 x tan x ^r ~ 33) degrees. In the top section, the 6 
cm PAs are as observed, not corrected for Faraday rotation. For comparison with the 2.2 /im PAs in the bottom section, 

a correction for a rotation of —106 rad m~ 2 has been applied. 
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Predicted 
2\xm Polarization 


Li 





6cm Polarization 



Fig. 1. — Total intensity and polarization for Cassiopeia A at 2.2 (im and 6 cm. The left side images of total intensity are at resolutions 
of 1" (approximate seeing) and 1.5", respectively. The center of all images is at RA 23h 23m 26s , Dec 58° 48' 58.4", (J2000), and the 
scale is shown on the 2.2 fim map. The images on the right arc of polarized flux, convolved to 10", along with the angle of the electric field 
(observed PA) for the boxes listed in Table 1. The middle frame on the right shows the predicted 2.2 fim polarization predicted from the 
radio measurements, as described in the text. 
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Box 1 original 
Ei Box 1 predicted 
De-reddened 
□ Box 1 
A Box 2 
O Box 3 
X Box 5 
• Box 6 



log Frequency (Hz) compressed scale 

FlG. 2. — Flux densities for five of the boxes shown in Table 1, as a function of frequency. The flux densities have been multiplied by a 
factor of iA 75 to accentuate the curvature. The errors in the radio are smaller than the symbol sizes. The errors at 2.2fim are dominated by 
uncertainties in the extinction correction. The open and closed squares and the dashed line connecting them show the size of the nominal 
extinction correction applied to each box. The error bar on the right shows the range of possible extinction corrections that could have been 
applied to the square (Box 1), as per the text discussion. The embedded squares show a prediction for Box 1 based on non-linear acceleration 
at modified shocks, again as described in the text. 




Azimuth (degrees) 



Fig. 3. — Inferred magnetic field angles (observed PA + 90 deg) and filament position angles as a function of azimuth. A characteristic 
error is shown for the 2.2 fim polarization and optical filament angles. The 6 cm angles have been corrected for the average RM in this region 
of -106 rad/m 2 and have errors of approximately the symbol size. The line indicates pure radial alignment. 



